Charge Density Wave Instability and Soft Phonon in APtsP (A^Ca, Sr, and La) 
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The electronic and phonon properties of the platinum pnictide superconductors APtsP (A=C&, Sr, 
and La) were studied using first-principles calculations. The spin-orbit coupling effect is significant 
in LaPtaP but negligible in CaPtaP and SrPtsP. Moreover, SrPtsP has been demonstrated to 
exhibit an unexpected weak charge- density- wave (CDW) instability which is neither simply related 
to the Fermi-surface nesting nor to the momentum-dependent electron-phonon coupling alone. The 
instability is absent in CaPtsP and can be quickly suppressed by the external pressure, accompanied 
with decreases in the phonon softening and BCS T c . Our results suggest SrPtsP as a rare example 
where superconductivity is enhanced by the CDW fluctuations. 
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I. INTRODUCTION 

The discovery of iron-based high-Tc superconduc- 
tors has stimulated considerable interests in the con- 
densed matter physics society. Sharing the similar 
layered square lattice crystal structures, both iron- 
pnictides/chalcogenides and cuprates exhibit the long- 
range magnetic order in parent compounds and the un- 
conventional superconductivity (SC) upon doping 1 6 . In 
these 3d transition-metal based materials, the electron 
correlation effect is considered to play a role in the for- 
mation of the spin-density- wave (SDW) instability while 
the electron pairing may be mediated or enhanced by the 
spin fluctuations. 

Recently, several families of platinum-pnictide based 
superconductors(SrPtAs 7 , SrPt 2 As 2 8 , APt 3 P 9 ) have 
been reported. Among these compounds, APtaP (A=Ca, 
Sr, and La) have shown considerable peculiarity. The 
crystal structures of these bd electron-based pnictides 
are experimentally determined to be anti-pevroskite with 
y/2 x y/2 x 1 distortion (FIG. 1). Due to the distortion, 
the platinum atoms take two different sites Pt 1 and Pt n , 
while the Pt 11 and P atoms form a Pt 2 P layer that re- 
sembles the iron-pnictide layers in the iron-based super- 
conductors (FIG. 1(b)). A particularly interesting fea- 
ture of this family is the significant variation of T c , i.e., 
T c = 1.5K, 6.6K, and 8.4K for A =La, Ca, and Sr, re- 
spectively. Among them, SrPtsP seems to be very spe- 
cial, as it has not only the highest T c in the 5d-electron 
based superconductors, but also a large 2Ao/ ksT c ~ 5, 
manifesting strong-coupling BCS superconductivity. It 
also exhibits a substantially lower electron specific heat 
coefficient 7 than its homologue CaPtaP and a notable 
anomaly near 10K slightly above T c . 

There are two major distinctions between the 3d and 
5d transition metal pnictides. First, the 5<i-orbitals are 
usually more extended so that the Coulomb interaction is 
relatively weak. Second, the spin-orbit coupling (SOC) 
in 5d-electrons is much stronger, as exemplified by the 
chemically similar compound SrPtAs 7 . In order to under- 
stand the puzzling features observed in the AS^P family, 
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FIG. 1: (color online) Crystal structure of APtsP compounds, 
(a) Crystal structure of APtsP. Yellow atoms are A (Ca, Sr, 
or La) atoms; red atoms are platinum; and blue atoms are 
phosphorus, (b) Pt 11 and P atoms form a Pt 2 P layer that 
is geometrically similar to the iron-arsenic layer in the iron- 
pnictides. Larger blue atoms denote the P atoms above the 
Pt plane and smaller blue atoms are the ones below the Pt 
plane. 



it is necessary and informative to clarify the spin-orbit 
coupling (SOC) and the electron-phonon coupling effects 
in the electronic structure of these materials. 

In this article, we present our latest first principles re- 
sults of APtsP (A=Ca, Sr, and La). We show that while 
spin-orbit coupling effect is negligible in SrPtsP and 
CaPtaP, it significantly affects the electronic structure 
of LaPtsP around the Fermi level Ef- An unexpected 
weak charge-density- wave (CDW) instability develops 
in the SrPtsP which is neither simply related to the 
Fermi-surface nesting nor to the momentum-dependent 
electron-phonon coupling alone. The instability is ab- 
sent in CaPtaP and can be quickly suppressed by the ex- 
ternal pressure, accompanied with decreases in the BCS 
T c . The CDW also reduces the electron density of states 
(DOS) around Ep for SrPtsP, leading to the anomalous 
specific heat behavior around 10K. 

The rest of this paper is organized as follows. In the 
next section, we briefly introduces the calculation method 
and details, followed by a section where we present in de- 
tail the electronic structure of APtsP (A=Sr, Ca, and 
La), the phonon properties of SrPtsP compared with 
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SrPt 3 P 


CaPt 3 P 


LaPt 3 P 


a (A) 


5.8788 (5.8852) 


5.7264 (5.7284) 


5.8187 (5.8221) 


c(A) 


5.4313 (5.4339) 


5.4507 (5.4618) 


5.5230 (5.5273) 


Z Pt 


0.1376 (0.1394) 


0.1370 (0.1382) 


0.1435 (0.1443) 


Z P 


0.2767 (0.2754) 


0.2686 (0.2672) 


0.2717 (0.2707) 



TABLE I: Optimized geometry of APtsP. z\> t denotes the frac- 
tional ^-coordinate of Pt 1 site. The numbers outside (inside) 
the brackets are obtained without (with) SOC. 



CaPt 3 P, as well as electron-phonon interactions and the- 
oretical T c of SrPt 3 P. We finally summarize and draw 
conclusions in section III. 



II. METHOD AND CALCULATION DETAILS 

To study the electronic structure of APt 3 P, we em- 
ployed the plane -wave basis projected augmented wave 
density functional (DFT) method implemented in the 
Vienna ab-initio simulation package (VASP) 10 ' 11 . The 
phonon properties were obtained using the frozen phonon 
method implemented in PHONOPY 12 . To ensure the 
calculation convergency to lmeV/cell, a 480 eV energy 
cutoff to the plane- wave basis and a 8 x 8 x 8 Monkhorst- 
Pack k-grid were employed to perform structural relax- 
ation until the maximum force on individual atoms is 
smaller than 0.001 eV/A and the internal stress smaller 
than 0.01 kbar. The DOS results are obtained using 
16 x 16 x 16 T-centered k-grid and the tetrahedra method. 
For electron-phonon coupling calculations and supercon- 
ducting T c estimation, we used density functional per- 
turbation theory implemented in the Quantum Espresso 
package 13 , where the electron-phonon coupling constants 
were evaluated with 4x4x4 T-centered phonon q-grid 
and 32 x 32 x 32 T- centered k-mesh. 



III. RESULTS AND DISCUSSION 

A. Geometry and Electronic Structures 

The optimized structural parameters are summarized 
in TABLE I, which agree very well with the experimen- 
tal results 9 (<3% difference). Thus the application of 
the DFT(+SOC) method is guaranteed in these mate- 
rials with a relatively weak Coulomb interaction. We 
then explore the electronic structures shown in FIG. 2. 
For both SrPt 3 P and CaPt 3 P, the SOC does not affect 
the electronic structure around the Fermi energy Ep, as 
no band splitting due to the SOC is observed within 
~ i^F±500 meV range. Hence the Fermi surfaces and 
consequently the SC of these two compounds are not af- 
fected by the SOC effect. The situation is, however, quite 
different in the case of LaPt 3 P, as it exhibits a signifi- 
cant band splitting due to the SOC effect from M to X 
around Ep, where an extra Fermi surface sheet emerges. 



This observation accounts for the significant reduction of 
T c in LaPt 3 P because the SOC-induced band splitting 
breaks the spin symmetry and thus reduces the pairing 
strength of the electrons. Consequently, LaPt 3 P may 
be a new member of the unconventional superconductors 
like SrPtAs 14 which deserves further experimental study. 
In the rest of our discussions, we will mainly focus on 
SrPt 3 P and CaPt 3 R 

Similarities between their electronic structures are ex- 
pected, since both Sr and Ca are iso-valent alkaline earth 
elements with two outmost s-electrons. Projected density 
of states (PDOS) plots (FIG. 2(d)-(e)) suggest that the 
DOS around Ep in both compounds are dominated by 
Pt-5d and P-2p orbitals, while the contribution from Sr 
or Ca orbitals are minimum. However, the band struc- 
tures of SrPt 3 P and CaPt 3 P show two significant differ- 
ences around Ep. First, one of the hole bands crossing 
Ep around the Z point in SrPt 3 P is missing in CaPt 3 P. 
Thus CaPt 3 P has one less Fermi surface sheets than 
SrPt 3 P(FIG. 3). Second, two bands cross Ep at approx- 
imately the same k point in SrPt 3 P from T to M, while 
these crossing k points are well separated in CaPt 3 P. 



B. Phonon Properties and Charge- Density- Wave 
Instability 

We further investigate the mechanical properties of 
these materials by performing phonon calculations us- 
ing the frozen phonon method. For SrPt 3 P, the re- 
sulting phonon band structure and DOS exhibit a clear 
instability evidenced by the appearance of the imagi- 
nary phonon frequencies. Remarkably, such instabil- 
ity is absent in CaPt 3 P. Usually, two scenarios can be 
responsible for this instability, namely, the SDW and 
CDW, respectively. The former appears frequently in 
the 3<i-transition-metal pnictide superconductors. We 
performed calculations with four possible spin configu- 
rations including/excluding the SOC. All these calcula- 
tions eventually converge to the non-magnetic spin con- 
figuration. Thus we exclude the SDW instability in this 
system. 

Noting that the soft phonon modes are located around 
the X (-7T, 7r, 0) and R (7r, 7r, tt) points (FIG. 5(a)), we 
then used a y/2 x y/2 x 1 supercell with respect to its orig- 
inal unit cell to adapt the lattice distortions due to the 
possible CDW instability. A less than 0.2% distortion to 
the atomic coordinates with y/2 x y/2 x 1 modulation is 
found to lower the total energy by ~ 10 meV/cell (or ~ 5 
meV/formula). It is also confirmed that the soft phonon 
mode disappears if the distorted supercell is used as the 
static initial condition to calculate phonon properties. 
Therefore, we demonstrated that the soft phonon modes 
are in fact due to the CDW instability in SrPt 3 P. It is 
important to emphasize that the energy scale of the pre- 
dicted CDW instability is about ~ 1 meV/atom, which 
is at least one order smaller than the SDW energy scale 
in iron pnictides. This indicates the CDW temperature 
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FIG. 2: Band structure and DOS of (a)(d) SrPt 3 P, (b)(e) CaPt 3 P, and (c)(f) LaPt 3 P at ambient pressure (0 GPa) with and 
without SOC. (a)-(c) are band structures, where red solid lines and blue dashed lines are obtained without and with SOC, 
respectively, (d)-(f) show DOS and its projection onto Pt-5d orbitals, with upper panels and lower panels showing the results 
with and without SOC, respectively. 




(a)Sr Sheet 1 (b)Sr Sheet 2 (c)SrPt 3 P 




(d)Ca Sheet 1 (e)Ca Sheet 2 (f)CaPt 3 P 

FIG. 3: (color online) Fermi surface sheets of (a-b) SrPt 3 P 
and (d-e) CaPt 3 P at GPa. Fermi surface cross section of 
(e) SrPt 3 P and (f) CaPt 3 P at k z = 0.0 (red) and k z = 0.5 
(blue). The centers of the squares are T(Z) and the corners 
are X. Notice that CaPt 3 P has one less fermi surface sheet 
around Z. 



is roughly around 10 K or even lower in SrPtsP. Fur- 
ther calculations show that the weak CDW instability is 
greatly suppressed by applying an external pressure and 




q-points 

FIG. 4: (color online) Fermi surface nesting function (/(q), 
red solid lines) and electron response function (xo(q), blue 
dashed lines) for SrPt 3 P at OGPa (top panel) and 5GPa (mid- 
dle panel), as well as CaPt 3 P at OGPa (bottom panel). 



disappears above 3 GPa. 

We next investigate the correlation between the CDW 
instability and the crystal structure by showing the evo- 
lution of the phonon DOS of SrPtsP under the external 
pressure(FIG. 5). At zero GPa, a DOS hump below 
THz is clearly presented due to the soft phonon modes 
caused by the CDW instability as shown above. The first 
real phonon DOS peak is formed around 2 THz, primar- 
ily contributed by the lowest optical modes. The exter- 
nal pressure compresses the interatomic distance, causing 
stiffer bonds, and therefore all real phonon modes shift to 
higher energies. For the soft modes, the DOS hump below 
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q-points 

(a)SrPt 3 P OGPa 




q-points 

(b)SrPt 3 P 3GPa 




Frequency (THz) 

(c)Phonon DOS 



FIG. 5: (a)-(b): Phonon band structure of SrPt 3 P at (a) 
GPa and (b) 3 GPa. (c): Phonon density of states under dif- 
ferent pressures. Negative frequencies indicate the imaginary 
(soft) phonon modes caused by CDW instability. 
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THz is reduced at 1 GPa and eventually disappears at 3 
GPa, indicating that the CDW instability can be quickly 
suppressed by the external pressure. The hardened soft 
modes contribute to the states below 2 THz, forming the 
new first peak below 2 THz. 

To analyze the formation of the CDW instability, we 
calculated the Fermi surface nesting function /(q) = 
Z)kmn^( e k+qn)^(ekm) and the non-interacting elec- 
tron response function X o(q) = ^ Ekmn for 
SrPtaP at different pressures, as well as for CaPtsP at 
GPa (FIG. 4) 26 . In all these cases, no significant di- 
vergence appears in both functions. Instead, /(q) shows 
peaks not only around X, but also around M, and Z. 
In the cases of SrPt 3 P at 5 GPa and CaPt 3 P at GPa, 
the nesting functions show even stronger peaks around 
Z. Besides, Xo(q) shows broad humps instead of sharp 
peaks around X, and that feature does not show much 
variation. Thus the Fermi surface nesting effect alone 
cannot be responsible for the CDW instability. Recently, 
it was suggested that in the dichalcogenides the CDW 
could be a result of strong electron-phonon coupling at 
certain q vectors 15 17 . If this is the case here, the q- 
resolved electron-phonon coupling constant A q of SrPtsP 
should exhibit an apparent anomaly around the CDW 
wave vector at GPa, which would then be suppressed by 
external pressure and disappear beyond 3 GPa. In order 
to clarify this possibility, we performed first-principles 
electron-phonon calculation at different pressures (FIG. 
6(a)). The resulting wave vector, however, is around 
(7r/2, 7r, 0), far from the predicted CDW vector (71", 7r, 
0). Therefore, the electron-phonon coupling alone can- 
not be accounted for the present CDW instability either. 



FIG. 6: (a) q-dependent electron-phonon coupling constant 
A q of SrPtsP under different pressures. (b)SC transition tem- 
perature T c of SrPtsP at different pressures. The circles in (b) 
indicates the T c obtained by first-principles electron-phonon 
coupling calculations and the Allen-Dynes formula; and the 
dashed line is a quadratic curve fitted to the results. 



C. Superconducting T c and Specific Heat Anomaly 

We now discuss the T c of SrPt 3 P within the BCS the- 
ory. A direct calculation of T c is prohibited due to the 
presence of the CDW-induced imaginary phonon modes. 
We thus performed a series of calculations at the pres- 
sures from 3 to 9 GPa, where the CDW instability is 
absent. The variation of the BCS T c under the ex- 
ternal pressure p is obtained by using the Allen-Dynes 
formula 18 , and then extrapolated to p = with least 
square fitting method assuming a quadratic dependence 
of T c onp(FIG. 6(b)). Hence the extrapolation represents 
the T c without the CDW instability. The resulting tran- 
sition temperature is 5.4 K, much lower than the 8.4 K 
as observed in experiment. The large difference between 
the calculated and experimental T c may be an evidence 
for CDW-enhanced superconductivity. 

In addition to T c , the experiment reveals that the 
specific heat coefficient 7, which is proportional to the 
electron DOS at the Fermi level g(Ep), is relatively 
larger in CaPtsP. Our calculation shows that in the ab- 
sence of the CDW, g(E F ) is 4.58 eV^/ceW for SrPt 3 P, 
which is slightly larger than the corresponding value 4.53 
eVVcell for CaPt 3 P. We then performed the DOS cal- 
culations using a y/2 x y/2 x 1 super cell that accom- 
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modates CDW distortions. The resulting g(E F ) = 4.31 
eV _1 /cell, which is 6% smaller than the one without 
CDW. Thus the formation of the CDW indeed leads to 
a sizable reduction in 7, in agreement with the exper- 
iment. The experimental data also reveal other delica- 
cies related to the specific heat and resistivity in SrPtaP 
slightly above the T c . While these features remain ques- 
tionable, we note that these features are just within the 
energy scale ~ 10K where the CDW instability starts to 
play a role. 

Several remarks are then in order. First, the CDW in- 
stability predicted for SrPtaP is rather unexpected. This 
is because of not only the relatively high T c ~ 8.4K ob- 
served in the 5d transition metal compounds, but also 
the very small CDW energy scale T* ~ 10K compared 
to the canonical CDW superconductors such as in a se- 
ries of transition metal dichalcogenides TiSe2 19 21 and 
NbSe2 22 ' 23 , as well as another platinum pnictide com- 
pound SrPt2As2 8 . So far the minimal difference between 
the CDW and SC temperatures in known CDW materi- 
als is observed in NbSe 2 with T c = 7.2K and T* = 32K. 
Second, while the CDW order competes and sometimes 
co-exists with the SC in most of known CDW supercon- 
ductors, the CDW in these cases is static, in the sense 
that the charge redistribution and ionic modulation are 
already stabilized and occur far above the SC instability. 
In the present case, however, the CDW instability should 
be considered as dynamical, because only a very small dis- 
tortion (less than 0.2%) in atomic coordinates is required 
to avoid the imaginary phonon modes and accommodate 
the lattice modulation, resulting in a very small charac- 
teristic temperature (~ 10K). Perhaps this is the reason 
why the CDW in SrPtsP has not been reported so far 
experimentally. It is also possible that quantum fluctu- 
ations may influence the weak CDW, resulting in quan- 
tum paraelectric states, similar to the ferroelectric state 
of SrTiOs. Third, in contrast to the conventional mech- 
anism of the CDW, neither pure nesting effect nor the 
momentum-dependence of the electron-phonon coupling 
can account to the CDW formation in SrPtsP indepen- 



dently. This implies that an appropriate combination of 
the nestings, their frustrations, the electron-phonon cou- 
pling, or even the states far away from the Fermi-surface, 
is at work in stabilizing the CDW in realistic materials 24 . 
We argue that this may be a characteristic feature of the 
dynamical CDW. 



IV. CONCLUSION 

In conclusion, we have performed a detailed study on 
APtaP using the first-principles simulations. Our calcu- 
lations show that the Fermi surfaces of these materials 
consist of multiple sheets. The spin-orbit coupling effect 
is negligible in both SrPtsP and CaPtaP but plays an im- 
portant role in LaPtaP. We predict that SrPtsP exhibits 
a CDW instability which can be quickly suppressed by 
the external pressure and is absent in CaPtaP under am- 
bient pressure. The formation of the CDW is neither 
directly related to the nesting effect nor to the electron- 
phonon coupling alone. Our results suggest that SrPtsP 
is a rare material where the SC instability falls into the 
fluctuation regime of the dynamical CDW which in turn 
leads to the enhancement in T c . 
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Note added. During the preparation of this 
manuscript, we became aware of another work by I. A. 
Nekrasov and M. V. Sadovskii 25 on the electronic struc- 
ture of SrPtsP. However, they did not consider the spin- 
orbit coupling effect, the phonon properties, nor the pos- 
sible CDW instability in this material. 
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